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A possibility to greatly enhance frequency-conversion efficiency of stimulated Raman scattering 
is shown by making use of extraordinary properties of three-wave mixing of ordinary and backward 
waves. Such processes are commonly attributed to negative-index plasmonic metamaterials. This 
work demonstrates the possibility to replace such metamaterials that are very challenging to engi- 
neer by readily available crystals which support elastic waves with contra-directed phase and group 
velocities. The main goal of this work is to investigate specific properties of indicated nonlinear 
optical process in short pulse regime and to show that it enables elimination of fundamental detri- 
mental effect of fast damping of optical phonons on the process concerned. Among the applications 
is the possibility of creation of a family of unique photonic devices such as unidirectional Raman 
amplifiers and femtosecond pulse shapers with greatly improved operational properties. 



I. INTRODUCTION 

Extraordinary features of coherent nonlinear optical 
NLO energy conversion processes in negative-index meta- 
materials (NIMs) that stem from wave-mixing of ordinary 
and backward electromagnetic waves (BEMW) and the 
possibilities to apply them for compensating optical losses 
have been investigated in |l|-[lO|. Essentially different 
properties of three-wave mixing (TWM) and four-wave 
mixing processes on one hand and second harmonic and 
third harmonic generation have been revealed in [l|, [lll - 
[la | . Ultimately, it was shown that NLO propagation pro- 
cesses that involve BW enable a great enhancement of 
energy-conversion rate at otherwise equal nonlinearities 
and intensities of input waves. A review can be found in 

In [1^ , three- wave mixing (TWM) of continuous waves 
(CW) was investigated in the scheme of stimulated Ra- 
man scattering (SRS) as an analog of parametric interac- 
tion of waves in a medium with negative dispersion. Two 
of the coupled waves were ordinary electromagnetic waves 
and the third one was backward elastic wave correspond- 
ing to optical phonons which exhibit negative dispersion 
duiv{k)/dk. The latter gives rise to contra-directed group 
and phase velocities of propagation of lattice vibrations, 
i.e. backward elastic waves. Such waves exist in crys- 
tals containing more than one atom per unit cell. Indeed, 
backwardness is a main property of EM waves propagat- 
ing in NIMs. Such BWs were predicted by L. I. Mandel- 
stam in 1945 [20], who also had pointed out that nega- 
tive refraction is a general property of the BWs. Optical 
phonons with wave vector much smaller than the recipro- 
cal lattice vector have been considered in [1^ . The focus 
was placed on a coupling scheme where the fundamen- 
tal radiation was converted to co-propagating Stokes and 



contra-propagating vibrational waves under the phase- 
matched conditions. The possibility of the extraordinary 
SRS was shown, which is typical for TWM whereas one 
of the coupled waves possesses negative dispersion [l|, [2[ ■ 
The latter results in gain of the Stokes component along 
the medium much greater than the exponential growth. 

The basic idea underlying the proposed concept is as 
follows. The dispersion curve uj{k) for optical phonons is 




FIG. 1. Negative dispersion of optical phonons and two phase 
matching options for short- and long-wave vibrations: (a) - co- 
propagating, (b) - contra-propagating fundamental, control, 
and Stokes, signal, waves. Insets: relative directions of the 
energy flows and the wave- vectors. 



depicted in Fig. [TJ It is negative in the range from zero 
to the boundary of the first Brillouin's zone. Hence, the 
group velocity of such phonons, vg*", is antiparallel with 
respect to its wave- vector, k^'', and phase velocity, v^'', 
because 



jU, 



(k/k)[duj/dk], duj{k)/dk < 0. (1) 
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Optical vibrations can be excited by the light waves due 
to the two-photon (Raman) scattering. The latter gives 
the ground to consider such a crystal as the analog of 
the medium with negative refractive index at the phonon 
frequency and to employ the processes of parametric in- 
teraction of three waves, two of which are ordinary EM 



waves and the third is the wave of elastic vibrations with 
the directions of the energy flow and of the wave- vector 
opposite to each other. For continuous wave SRS in the 
coupling geometry depicted in Fig. [IJa) and diamond 
crystal, estimates made in tlSl] have shown that the min- 
imum intensity of the fundamental field /min, required to 
reach the aforementioned extraordinary Raman amplifica- 
tion is on the order of /mm "^ 10^® W/cm^, which exceeds 
optical breakdown threshold. The main process that de- 
termines such a high intensity is fast phonon damping, 
which is characterized by the optical phonon relaxation 
rate on the order of r„ = 6 • 10^^^ s, and relatively small 
group velocity of the elastic oscillation Vy ^ 10^ -^ 10'^ 
cm/s for optical phonons with small optical wave vectors 
|2ll . 1221 . Diamond is a representative of a wide class of 
transparent crystals which support optical phonons [21- 

The goal of this work is to show that this seemingly 
formidable obstacle can be removed and the threshold in- 
tensity of the fundamental field /min required to realize 
such extraordinary coupling can be significantly reduced. 
It is possible by making use of short laser pulses with 
duration Tp less then the lifetime of the vibrational oscil- 
lation of Ty. Some of the properties of the output fun- 
damental and Stokes pulses, such as duration and pulse 
shape are investigated. 



II. EQUATIONS AND MODEL 

Consider lowest-order SRS process [2^, [23l ■ The waves 
are given by equations 

S,., = (l/2)e,,,(z,t)e*-^~'"'-* -fee, 
Qy = (l/2)Q(z, t)e*'="--'""* + c.c. (2) 

Here, ei^s, Q, ^i,s.v and ki,s,y are the amplitudes, frequen- 
cies and wave-vectors of the fundamental, Stokes and vi- 
brational waves; Qy(z,t) — ,ypx{z,t); x is displacement 
of the vibrating particles, p is the medium density. With 
account for the frequency and phase matching exppres- 
sions, 

UJl = UJs + LOy (ky) , h^ks [Uls] + ky, 

one obtains following partial differential equations for the 
slowly varying amplitudes in the approximation of the of 
first order of Q in the polarization expansion 28;] : 

dz vi dt kic^ dQ 
d£s 1 d£s .TTwJ ,,9a 



az Vg at kgC' oQ 
dQ ^ j_dQ ^ Q ^ , 1 .rdc^c^c* 

dz Vy dt ly 4:U!yVy OQ 



(3) 



Here, vi^s.v are the projections of the group velocities of 
the fundamental, Stokes and vibration waves on the z- 
axis, N is the number density of the vibrating molecules, 



a is the molecule polarizability, ly — TyVy is the mean free 
path of phonons. 

For the sake of simplicity, we consider model of a rect- 
angular pulse of input fundamental radiation with dura- 
tion much shorter that the phonon lifetime Ty. In the 
coordinate frame associated with this pulse and within 
the fundamental pulse range, complex amplitudes of two 
other interacting fields become time independent, and 
Eqs. ([3]) transform into set of ordinary differential equa- 
tions whose solution is known in the approximation of 
constant pump amplitude or can be relatively easily found 
numerically when depletion of the pump is accounted for. 
Here, the boundary conditions must be fulfilled not at the 
boundaries of the medium but at the boundaries of the 
fundamental pulse. The latter is correct for the period 
of time after the instant when generated waves reach the 
boundaries of the fundamental pulse due to the differ- 
ence in their group velocities or direction of propagation. 
Such approximation becomes true after travailing a dis- 
tance I > L^'^"", where L'^'''' = max{L^,L^}. Here, L^-", 
which is further referred to as group length, is defined as 



'Tpvf/lVs 



Vl\ 



(4) 



where Vg^y are negative if their energy fluxes are oppo- 
site to that of the fundamental wave. Hereinafter, the 
waves are referred to as co-propagating waves if Poynt- 
ing vector of the Stokes wave is co-directed with that of 
fundamental wave (Fig. [T^), and as counter- propagating 
in the opposite case (Fig. [T})). In the approximation of 
constant pump amplitude, in the coordinate frame locked 
to the pump pulse and within the pulse, equations for 
generated Stokes and backward vibration waves take the 
form: 



dQ/dS, = -igv£*s - QKy/ly 

dSs/dS, = igsQ*, 



(5) 



where gy = -KyN{da/dQ)£i/{ALUyVy), Ky^g = 
Vv.s/{vy,s - vi), gs = KgN[da/dQ)£i-Kijj1/{ksc'^). Here, 
intensity of the laser beam is taken constant. Depletion 
of the fundamental beam is neglected. Equations ([5]) 
may describe significant amplification of the Stokes sig- 
nal, however they remain valid until only relatively small 
part of the strong input laser beam is converted. Equa- 
tions ([S]) do not depend on time and are similar to those 
describing CW TWM 1^. Here, group velocities Vy < 0, 
i;^ > and wave vector kg > for co-propagating and 
Vg < 0, kg < for counter-propagating phase matching 
configurations. Since the Stokes frequency is less than 
that of the fundamental one, Vg > vi and is much greater 
than the magnitude of phonon group velocity Vy.ln the 
case of co-directed laser and Stokes waves, the boundary 
conditions take the form: 



f,(e-o)-£°, Q(e = g-o. 



(6) 



In the opposite case of counter-propagating Stokes and 
laser waves, the boundary conditions are given by equa- 
tions 



£gi^^ lp)^£l^, Q{^^lp)^0. 



(7) 



Here, Ip — TpVi is length of the fundamental pulse. 
Solution to Eqs. ([5]) can be presented in the form: 



S 



p7'?' 



{Ri cos Z + y sin Z} 



Ri cos (Ri) + y sin (Ri) 



2R2 



(8) 



Here, T^^ = |f.(0/f°|' and rti = |f.(e)/d"| are 
transparency (amplification) factors, 7' = 7^^, 7 = 

/3i,2 = 7 ± i?2, .g' = g/p, g = ^.g^^ffs- 

For 7 <C g, phonon damping can be neglected and 
Eqs. © take the form 






cos [g?p(l - C)] 



cos (^Zp) 

'p(i-e') I 



,-glpil-e 



/2 



(9) 



III. ENHANCING COHERENT ENERGY 

TRANSFER BETWEEN ELECTROMAGNETIC 

WAVES THROUGH BACKWARD OPTICAL 

PHONONS 

In the given approximation of neglected depletion of 
the fundamental wave, solution to Eq. (|9]) tends to infin- 
ity for certain pulse length and intensity. This indicates 
the possibility to greatly enhance the conversion efficiency 
provided that the following requirements are met: 



'S (glp) =0; lp = {TT/2)/y/^, 



(10) 



Corresponding threshold intensity of the fundamental 
field /^in, which is required for realization of great en- 
hancement of energy conversion due to NLO coupling 
with BW phonon, is given by equations: 



5 > 7, ^mi 



Kg IQlT^VyT^ 



N 



da 
dQ 



(11) 



Here, Hs is refractive index at Stokes frequency, Aso is 
Stokes wavelength in the vacuum. 



A. Factors discriminating BW SRS in continuous 
wave and pulsed regimes 

Threshold intensity of fundamental radiation required 
for realization of BW SRS in CW regime is given by JK| 

/min = {cnAsoi^v/Sn%T) \Nda/dQ\-^ . (12) 

It differs from the threshold value in pulse regime by fac- 
tor 



I" K 



Vy Vs - Vl 



(13) 



For the same typical crystal parameters as employed in 
0, Eq. dni) yields /^,i„//mi„ « 10-". Hence, I^,^ de- 
creases down to /^jj^ ~ 10^ W/cm^, which is achievable 
with commercial femtosecond lasers and falls below opti- 
cal breakdown threshold for most transparent crystals. 

Equation (jl3|) displays two factors that determine sub- 
stantial decrease of /^j„ in pulsed regime compared to 
that in CW one. First factor is the ratio of group velocity 
of the elastic wave to that of the fundamental one, Vy/vi, 
which is on the order of ~ 10"^. This factor is attributed 
to the fact that phonons generated on the front edge of 
the laser pulse propagate in the opposite direction and, 
hence, exit very fast, practically with the optical group 
velocity, from the fundamental pulse zone before it is dis- 
sipated. Hence, effective phonon mean free pass becomes 
commensurable with the fundamental pulse length. This 
mitigates the detrimental effect of phonon damping. With 
increase of Vy, phonon mean free path grows, which de- 
creases both /^jj^ and /min in a way that the advantage 
of pulse regime over CW regime diminishes. The second 
factor in Eq. ([13]) determines further decrease of I^^^ due 
to small optical dispersion in the transparency region of 
the crystals. The fact that the Stokes pulse surpasses the 
fundamental one slowly increases significantly the effec- 
tive NLO coupling length. 



B. Numerical simulations, Tp <g r„. 

Numerical analysis is done for the model with param- 
eters typical for diamond crystal [2l|, |22{ : carrying wave- 
length of the fundamental pulse A = 800 nm, pulse du- 
ration Tp = 60 fs, ujy = 1332 cm^^, vibrational transi- 
tion width {cTyy^ = 1.56 cm^i, vi = 1.228 • 10^° cm/s, 
Vs = 1.234 • 10^" cm/s, Vy — 100 cm/s for co-propagating 
and Vy = 2000 cm/s for counter-propagating waves, 
Nda/dQ = 3.78 • 10^ (g/cm)i/2. Partial differential equa- 
tions ([3]) were solved numerically in three steps: TWM in 
the vicinity of the entrance to the Raman medium, TWM 
and propagation through the medium and TWM in the 
vicinity of the exit from the Raman slab. Numerical sim- 
ulation for the first and third medium intervals were made 
in the laboratory reference frame with the boundary con- 
ditions applied to the corresponding edges of the slab. 
NLO propagation process inside the slab was simulated 
in the moving frame of reference with the boundary con- 
ditions applied to the pulse edges. Such a computing ap- 
proach allows for significant reduction of the computation 
time since, for each given instant, integration is required 
only through a space interval covered by the fundamen- 
tal pulse and not through the entire medium. Shape of 
the fundamental pulse was chosen nearly rectangular and 
symmetric with respect to its center 



£1 ^ -£,"{tanh[(to + tp 



t)/tf]~tanh[{to-t)/tf]}, 



the slope tf = 0.1, the pulse duration at half- maximum 
ip = 1 and pulse delay Iq = 0.6 were scaled to the fun- 
damental pulse width Tp. The amplitude of input CW 
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FIG. 2. Quantum conversion efficiency vs. intensity (energy) 
of tfie input pulses for co-propagating and contra-propagating 
pliase matcliing geometries. 



Stokes signal was chosen f ° = 10 ^£^ ■ 

Figure [5] displays output quantum conversion efficiency 
r\sq = {lvi/uJs) ■ J^Isiz,t)dt/ J^Ii{z = 0,t)dt vs. input 
pulse intensity (the same as vs. input pulse energy) 
both for co-propagating (z = L) and counter-propagating 
(z = 0) geometries. A great increase of the conversion ef- 
ficiency due to BW effect in the case of co-propagating 
waves is expficitly seen. Saturation at l)* / I^i^ ~ 7 • 10'' 
is due to depletion of fundamental radiation caused by 
conversion to Stokes radiation. 

Simulations also show that shape of the output Stokes 
and fundamental pulses differ and vary significantly de- 
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FIG. 3. Pulse sliape of amplified Stokes radiation at tlie 
exit of crystal of length L = 1 cm, Ts — \£s{L,t)/£s\ 
(solid line). Pulse shape of the output fundamental pulse 
Ti = jf,(L,t)/£:fj^isdepictedbydashedline. Here, jf/JP.„ = 
4 X 10*. 
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FIG. 4. Pulse shape of amplified Stokes radiation (solid line) 
and of transmitted fundamental radiation (dashed line) at the 
exit of crystal of L = 1 cm. ry^ — {loi/ljs) ■ \£s{L,t)/£J^\ is 
conversion efficiency, Ti = \£i{L,t)/£^\ is transmittance. a. 
I?/lLn = 8.3 X 10*. b. /f//L„ = 1.67 X 10^ 



pending on the intensity of the input fundamental wave. 
Input seeding signal at Stokes frequency is assumed a 
weak CW broadband radiation. Figure |3] depicts am- 
plified output pulse Ts = \£s{L,t)/£s\ of Stokes radi- 
ation for relatively small amplification of the signal and 
small depletion of the fundamental beam. Here, shape of 
the fundamental pulse remains unchanged. Pulse shape 
of amplified Stokes signal is different and determined by 
the fact that Vs > vi and Stokes pulse surpasses fun- 
damental one. In contrast, at I°/I^.in — ^-^ ^ 10^ ^^^d 
If/Imin — 1-67x10^ (Figs. 2^, b), conversion efficiency be- 
comes significant {rjq — 0.8 and rjq — 0.93, respectively). 
Corresponding depletion of the output fundamental pulse 
and changes in its shape is explicitly seen. Note that in 
the case shown in Fig.Ufb), the output Stokes pulse signif- 
icantly overtakes the pump pulse. In the latter case, ma- 
jor conversion occurs in the middle of the medium, after 
which both pulses propagate almost without interacting 



and independently from each other. It is seen that output 
Stokes pulse narrows with increase of the input funda- 
mental intensity. Here, crystal length of 1cm corresponds 
to more than 1000 input pulse lengths {T/tp = L/Lp — 
1357). Threshold intensity /^.„ = 6 • 10^ W/cm^, which 
corresponds to 60 fs pulse of 5 ^J focused to the spot of 
diameter D — 100 ^m. Intensity of seeding Stokes signal 
was chosen Pjl^ = 10"!°. 

Described NLO propagation processes are in striking 
contrast with their counterparts in crystals where only 
phonons with positive group velocity exist j2^, [23] . Such 
NLO properties are also different from those inherent to 
phase-matched mixing of EM and acoustic waves for the 
cases where the latter have energy flux and wave vec- 
tor directed against EM waves [20|. Elaboration of the 
proposed concept allows to utilize revealed extraordinary 
features for creation of a family of unique photonic de- 
vices made of ordinary Raman crystals such as optical 
switches, filters, amplifiers and cavity-free optical para- 
metric oscillators. Proposed here concept is different from 
earlier proposed in 30] and does not require periodic pol- 
ing of quadratic nonlinear susceptibility of crystals at the 
nanoscale as described in [31[ (and references therein). 



sponse by such mesoatoms also presents a challenging goal 
not yet achieved. This paper proposes to mimic similar 
extraordinary backward-wave extraordinary nonlinear- 
optical propagation processes, however, making use read- 
ily available Raman active crystals. Among the prospec- 
tive application is a family of photonic devices with ad- 
vanced functional properties. Basic underpinning idea is 
to replace one of the coupled backward electromagnetic 
wave by optical phonons - elastic waves with negative 
group velocity. Operation in short pulse regime is proposed 
to remove such a severe detrimental factor as fast phonon 
damping. Significant decrease of the required minimum 
intensity of fundamental radiation compared with that 
in continuous-wave regime is shown down to that pro- 
vided by commercial lasers. Unusual properties of the 
investigated processes are numerically simulated and the 
possibilities of pulse shape tailoring are predicted. There 
are many Raman active transparent crystals that support 
optical phonons with negative group velocity can be uti- 
lized, some of them, such as calcite, may offer further op- 
timization of operational characteristics of the proposed 
photonic devices. 
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